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Abstract - The Moderate Resolution Imaging Spectroradiometer (MODIS) is 
an earth- vie wing sensor that flies on the Earth Observing System (EOS) Terra 
and Aqua satellites, launched in 1999 and 2002, respectively. MODIS scans a 
swath width of 2330 km that is sufficiently wide to provide nearly complete 
global coverage every two days from a polar-orbiting, sun-synchronous, plat- 
form at an altitude of 705 km. MODIS provides images in 36 spectral bands be- 
tween 0.415 and 14.235 jum with spatial resolutions of 250 m (2 bands), 500 m (5 
bands) and 1000 m (29 bands). These bands have been carefully selected to en- 
able advanced studies of land, ocean, and atmospheric properties. Twenty-six 
bands are used to derive atmospheric properties such as cloud mask, atmos- 
pheric profiles, aerosol properties, total precipitable water, and cloud properties. 
In this paper we describe each of these atmospheric data products, including 
characteristics of each of these products such as file size, spatial resolution used 
in producing the product, and data availability. 


I. Introduction 


The Moderate Resolution Imaging Spectroradiometer (MODIS) is an EOS fa- 
cility instrument that is currently flying aboard the Terra and Aqua spacecraft. It 
is especially well suited to global monitoring of atmospheric properties from 
space, and is based on heritage sensors such as the Advanced Very High Resolu- 
tion Radiometer (AVHRR), Landsat Thematic Mapper (TM), High-resolution In- 
frared Radiation Sounder (HIRS), and the Nimbus-7 Coastal Zone Color Scanner 
(CZCS). The wide spectral range (0.41-14.24 fim), frequent global coverage (one 
to two days revisit), and two high spatial resolution bands (250 m), permit state 
of the art global monitoring of atmospheric profiles, column water vapor 
amount, aerosol particles, and the subsequently formed clouds [10]. Barnes et al. 
[2] provide a detailed description of MODIS, including its performance attrib- 
utes, optical design, spectral band characteristics, primary purpose of each band, 
pixel size, and signal-to-noise ratios at specified radiance levels, information that 
is essential for an in-depth understanding of the onboard calibrators and the op- 
eration of this highly sophisticated sensor. 

The intent of this paper is to provide brief descriptions of each of the MODIS 
atmosphere data products that are being produced by the Terra mission and to 
describe differences, where they occur, with comparable products to be pro- 
duced by Aqua. MODIS data are processed into various processing levels from 
Level- 1 (radiances and brightness temperatures that have been geolocated), to 
Level-2 (derived geophysical data products at the same resolution and location 
as the Level-1 data) to Level-3 (variables mapped onto uniform space-time grid 
scales). In this paper we describe the essential elements of all Level-2 'pixel- 
based' MODIS atmosphere data products, and give examples of selected parame- 
ters in each product. Finally, we describe the MODIS Level-3 atmospheric prod- 
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ucts that are produced at daily, eight-day, and monthly time intervals on a global 
1° x 1° latitude-longitude grid. In addition to simple statistics (mean, standard 
deviation, etc.) computed for each parameter, the Level-3 products also contain 
marginal density and joint probability density functions between selected pa- 
rameters. 


II. The MODIS Atmosphere Data Processing System 

Terra and Aqua data are stored on the spacecraft in solid-state memory and 
periodically transmitted to the ground. Once on the ground, these data are sepa- 
rated into instrument packets 2 hours in length and sent to the Goddard Space 
Flight Center for further processing. Fig. 1 shows the data processing architec- 
ture and data products that are used in processing MODIS atmospheric proper- 
ties (adapted from [20]). The first step is to produce Level-1 A (MODOl) data that 
consist of counts for all 36 bands of MODIS, along with raw instrument engi- 
neering and spacecraft ephemeris data. This data product is separated into 5 
minute granules, and hence all Level-IB, geolocated, and subsequently processed 
pixel-level data are processed into 5 minute granules (~2000 km in length along 
the orbital track), producing some 244 granules (or files) per day per data prod- 
uct. 

Following Level-IA, calibration is applied using the onboard calibrators, 
space looks, and instrument performance data, and these Level-IB (MOD02) data 
are stored in 3 separate files, one each for 250, 500, and 1000 m bands. Since 
many algorithms that follow in the subsequent processing chain require a com- 
bination of spectral bands having inherently different spatial resolutions, a Level- 
IB product is produced at 250 m resolution (daytime, 2 bands), 500 m (daytime, 5 
bands plus two 250 m bands aggregated to equivalent 500 m bands), and 1000 m 
(daytime, 36 bands all at 1000 m or aggregated to that resolution from the finer 
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resolution bands). At night, only the thermal bands are recorded so the MOD02 
1 km granules are smaller in size (17 bands) than during the daytime when all 36 
bands are recorded. Geolocation of each pixel is also contained in the MOD02 
file, but only every 5 th pixel and every 5 th scan line. For full resolution naviga- 
tion, a separate geolocation (MOD03) data file is produced at a 1 km spatial 
resolution, containing geodetic coordinates, ground elevation, and viewing ge- 
ometry. 

All of the data products described in this paper are Level-2 and Level-3 
products that ultimately require well calibrated and geolocated data from 
MOD02 and MOD03. These data products include the cloud mask (MOD35) for 
distinguishing clear sky from clouds, atmospheric profiles (MOD07), aerosol 
product (MOD04), precipitable water product (MOD05), cloud product 
(MOD06), and gridded time-averaged (Level-3) atmosphere products on a daily, 
eight-day, and monthly mean basis (MOD08 products). Table I lists the MODIS 
atmosphere data products, responsible investigators, processing level, product 
ID used in ordering data from the EOS data system, spatial resolution, file fre- 
quency, and file size. All of these data products are archived and distributed by 
the Goddard Distributed Active Archive Center (DAAC). Data products from 
Terra use the product ID MODxx, whereas those from Aqua use the product ID 
MYDxx (e.g., MYD06). 


III. Cloud Mask 

The cloud mask (MOD35) algorithm classifies each pixel as either confident 
clear, probably clear, probably cloudy, or cloudy. The cloud mask consists of 48 bits 
of output that include information on individual cloud test results, the process- 
ing path, ancillary information (e.g., land /sea tag), and additional information 
(e.g., thin cirrus detected, heavy aerosol). The cloud mask algorithm [1] uses a 


IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, OCTOBER 2002 


4 


series of threshold tests applied to 17 of the 36 MODIS bands to identify the pres- 
ence of clouds in the instrument field-of-view. The specific tests executed are a 
function of surface type, including land, water, snow /ice, desert, and coastal, 
and are different during the day and night. Each cloud detection test returns a 
confidence level that the pixel is clear ranging in value from 1 (high confidence 
clear) to 0 (low confidence clear). Tests capable of detecting similar cloud condi- 
tions are grouped together and a minimum confidence is determined for each 
group. The final cloud mask is then determined from the product of the results 
from each group. This approach is clear-sky conservative in the sense that if any 
test is highly confident that the scene is cloudy, the final clear sky confidence is 0. 
The first 8 bits of the cloud mask provide a summary adequate for many proc- 
essing applications. 

Fig. 2 shows an example of the cloud mask result for a daytime granule over 
Western Australia on June 5, 2001. The 'true color' image is a composite image of 
MODIS bands 1, 4 and 3 (0.645, 0.555, and 0.469 jim, respectively). The accom- 
panying cloud mask image shows the final result of the cloud mask algorithm 
with confident clear (blue), probably clear (green), probably cloudy (yellow), and 
cloudy (white). Most of the region is flagged as confident clear or probably clear. 
Some of the land near the coast is flagged as probably clear while the true-color 
composite suggests the region to be clear. This results from a combination of 
thin cirrus, haze, and varying surface radiative properties. Experience gained in 
comparing the cloud mask results with collocated surface lidar systems suggest 
that the confident clear and probably clear designations are invariably clear, 
whereas the probably cloudy and cloudy designations are invariably cloudy, 
with the probably cloudy designation arising near the edge of clouds or where 
there is thin cirrus present in the atmosphere. Subsequent algorithms are given 
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the results of the individual tests and may choose to use these to further deter- 
mine whether to process the pixel, as different algorithms show different sensi- 
tivities to these cloud mask results. 

IV. Atmospheric Profiles 

The atmospheric profiles (MOD07) algorithm performs statistical retrievals 
of atmospheric temperature and moisture layers, total precipitable water, total 
column ozone, and stability indices. While MODIS is not per se a sounding in- 
strument, the information content of the high spatial resolution infrared multi- 
spectral radiance observations can improve upon a priori definitions of atmos- 
pheric state by providing better delineation of horizontal gradients. The MODIS 
clear-sky retrievals are performed over land and ocean for both day and night 
when at least 20% of the radiances measured within a 5 x 5 field of view area 
(approximately 5 km resolution) are cloud-free. The retrieval methods employed 
here are partly based on the work of the International ATOVS (Advanced TIROS 
Operational Vertical Sounder) Processing Package (IAPP) [15] and the Geosta- 
tionary Operational Environmental Satellite (GOES) sounder algorithms [16], 
[13], [14]. 

The operational MODIS atmospheric profile algorithm is a statistical regres- 
sion. The algorithm uses 12 infrared bands with wavelengths between 4.47 pm 
(band 24) and 14.24 pm (band 36). Surface emissivity effects in the shortwave in- 
frared bands are mitigated by regressing against band differences (e.g., instead of 
T 4.52 and T 4.47 we use the difference T 4 . 52 -T 4.47 in the regression, where T>. rep- 
resents brightness temperature at wavelength k). An extension of the NOAA -88 
data set containing more than 8400 global radiosonde profiles of temperature, 
moisture, and ozone is used to determine the regression coefficients. MODIS in- 
frared band radiances are calculated from the radiosonde observations of the at- 
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mospheric state, and the regression coefficients are generated from these calcu- 
lated radiances/ atmospheric profile pairs. The radiative transfer calculation of 
the MODIS spectral band radiances is performed using a transmittance model 
called the Pressure layer Fast Algorithm for Atmospheric Transmittances [4]; this 
model has 101 pressure level vertical coordinates from 0.05 to 1100 hPa. The cal- 
culations take into account the satellite zenith angle, absorption by well-mixed 
gases (including nitrogen, oxygen, and carbon dioxide), water vapor (including 
the water vapor continuum), methane, carbon monoxide, nitrogen dioxide, and 
ozone. In addition, MODIS instrument noise is added into the calculated spec- 
tral band radiances. 

To perform the retrieval, the regression coefficients are applied to actual 
MODIS measurements to obtain the estimated atmospheric profiles, the integra- 
tion of which yields the total precipitable water (TPW) 1 or total column ozone. 
The advantage of this approach is that it does not need MODIS radiances collo- 
cated in time and space with atmospheric profile data, as it requires only histori- 
cal profile observations. To address radiative transfer model uncertainties and 
residual instrument calibration anomalies, radiance bias adjustments are made in 
the retrieval algorithm; these are calculated seasonally from MODIS observed 
clear sky radiance composites and global model atmospheric analyses. 

Fig. 3 shows the temperatures retrieved from the Terra-MODIS radiances on 
May 22, 2001 at 0500 UTC at levels (a) 850 hPa. (b) 700 hPa, and (c) 500 hPa over 
the Midwestern United States. At this time, the cold front associated with a low- 
pressure system centered in northern Wisconsin extended southward through 


1 Atmospheric profile retrievals are saved at only 20 levels in the MOD07 data product; however, 
the integration to TPW is performed from the original retrieved profiles with 101 levels. Thus, 
integration of the 20-level profiles will not necessarily result in the same value reported in the 
TPW field. 
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Illinois, Missouri, and eastern Texas. The tongue of cool air behind this front can 
be seen pushing into the central and upper Midwest, and the pattern extends up 
through all three levels in the atmosphere shown. 

Preliminary evaluation of the Aqua-MODIS atmospheric profiles algorithm 
(MYD07) indicates its performance is comparable to that of Terra. Fig. 4 shows 
one example from the Aqua infrared temperature and moisture retrievals. The 
total precipitable water (in cm) from September 4, 2002 at 1750 UTC indicates a 
sharp moisture boundary along the east coast of the United States associated 
with the cold front of a low-pressure system in eastern Canada. The associated 
GOES-8 Sounder precipitable water shows general agreement, but it does not 
have the spatial resolution to capture the sharp moisture boundary. A cross- 
section drawn through the strongest moisture gradient highlights the vertical 
structure of this front. Moist air, with mixing ratios up to 15 g kg' 1 , exists 
through 700 hPa to the south of the front, then the air dries out sharply around 
38.4° latitude. The vertical structure in the temperature field shows less structure 
than the moisture, and the temperature difference across the front is also less 
pronounced; however, warmer surface temperatures are evident at latitudes 
south of the front. 

Comparisons of MODIS atmospheric products with other observations collo- 
cated in time and space are facilitated by specialized instrumentation at the 
Southern Great Plains (SGP) Atmospheric Radiation Measurement-Cloud and 
Radiation Testbed (ARM-CART) site in Oklahoma. Terra and Aqua each pass 
over the SGP CART site twice daily; radiosondes are launched three times each 
day; and observations of total column moisture are made by the microwave ra- 
diometer (MWR) every 40-60 sec. Additional comparisons are possible with the 
GOES-8 sounder [18], [19] that retrieves TPW hourly. Comparison for 17 months 
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of Terra MODIS infrared TPW for 80 dear-sky cases from April 1 to September 1, 
2002 with TPW from GOES-8 sounder, radiosondes, and MWR are presented in 
Fig. 5. These cases were chosen by manual inspection to exclude any scenes with 
cloud contamination and with MODIS sensor zenith angle 0 > 50° at the CART 
site. MODIS shows general agreement with all three systems, the RMS differ- 
ence between MODIS and MWR TPW is 4.7 mm (compared with 1.4 mm and 1.3 
mm for GOES-8 and radiosondes, respectively). For dry atmospheres, MODIS 
overestimates the total column moisture; the average TPW bias (MODIS minus 
MWR) is approximately 3 mm for the 36 dry cases with TPW s 17 mm. On aver- 
age for the 44 wet cases (TPW > 17 mm), the MODIS regression-based TPW re- 
trievals are drier than the MWR by 5.7 mm. Thus MODIS will reduce the magni- 
tude of a TPW gradient; nonetheless, MODIS will still capture the small-scale 
structure of spatial gradients across a front as well as the global changes from 
day to day. 

MODIS is producing global fields of temperature, moisture, and ozone. The 
primary value of these fields lies in the depiction of fine scale horizontal gradi- 
ents. When combined with the increased atmospheric profiling capability of 
high spectral resolution infrared systems (such as the Atmospheric Infrared 
Sounder, AIRS), the MODIS offers improved sounding performance in cloudy 
conditions (by identifying sub AIRS field-of-view cloudiness) and delineates spa- 
tial gradients beyond the AIRS instrument resolution. 

V. Aerosol Product 

The aerosol product (MOD04) is based on different algorithms for the remote 
sensing of tropospheric aerosol over land [9] and ocean [25]. Both algorithms try 
to match MODIS observed reflectances to a lookup table of pre-computed reflec- 
tances for a wide variety of commonly observed aerosol conditions, as summa- 
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rized by King et al. [11]. Over land, the prime difficulty is separating the reflec- 
tance measured by the satellite into an atmospheric part and a land surface part. 
This difficulty is overcome by estimating the visible surface reflectance from the 
MODIS-measured reflectance in the shortwave infrared wavelength (2.13 gm) for 
dark targets, as described by Kaufman et al. [9]. Removing the estimated surface 
reflectance from the total measured by the satellite, the atmospheric contribution 
is isolated by using lookup tables to determine the aerosol characteristics and op- 
tical thickness in the scene. Final land products include aerosol optical thickness 
(t 3 ) at 0.47, 0.56 and 0.65 gm at a 10 km spatial resolution. The spectral depend- 
ence of the reflectance across the visible wavelengths is then used to obtain a 
rough estimate of the fine mode (radius < 0.6 gm) fraction of the aerosol optical 
thickness at 0.56 gm. 

The ocean algorithm assumes that the surface reflectance contribution to the 
total reflectance is small and can be calculated from geometry and assumptions 
about the sea state [25]. The better-characterized ocean surface permits the use of 
reflectances at 6 wavelengths (0.56-2.13 gm) in the derivation. Again, the re- 
trieved aerosol products are represented by the best fits between observed reflec- 
tance and the lookup table. Ocean products include aerosol optical thickness at 
0.47, 0.56, 0.65, 0.86, 1.24, 1.64, and 2.13 gm at a 10 km spatial resolution. Using a 
wide spectral range over ocean allows the derivation of the effective radius of the 
particle population, and more accurate determination of fine mode fraction than 
over land. The aerosol algorithm is applied to cloud free conditions outside the 
sunglint region (glint angle < 40°) and for land surface reflectance at 2.13 gm < 
0.25 at nadir and 0.4 at slant view directions. The comparison with the lookup 
table is performed for the average of the 25%-75% pixels over the ocean and 20%- 
50% over the land in order to avoid any undetected non-aerosol contamination. 



IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, OCTOBER 2002 


10 


A minimum of 10 (12) pixels is required for the analysis in any 10 km region over 
ocean (land). 

Fig. 6 shows a multispectral "true-color" composite image, associated re- 
trieved T a (0.56 jum), fraction of fine mode aerosol to total optical thickness (tf/x a ) 
at 0.56 jum, and aerosol effective radius for one granule at 16:25 UTC on May 4, 
2001. This granule, located over the east coast of the United States, displays con- 
sistency between the entirely independent land and ocean retrievals. Areas 
where the aerosol optical thickness cannot be retrieved due to the presence of 
clouds or sunglint are colored gray in Fig. 6b. Aerosol effective radius, shown in 
Fig. 6d, is only derived over oceans. The retrieval in this granule shows the aero- 
sol over land and near-shore coastal region to be dominated by fine mode aero- 
sols with an effective radius ~0.2 jum, typical of pollution episodes in this region 
[22]. Further from land the pollution gradually gives way to larger aerosol, pos- 
sibly a mixture with sea salt. The accuracy of satellite estimates of aerosol optical 
thickness was first suggested based on theoretical analyses [9], [25], and con- 
sisted of a bias of ±0.05 over land and ±0.03 over ocean due to uncertainty in the 
estimate of surface reflectance, and ±0.20t a over land and ±0.05t a over ocean due 
to uncertainty in the aerosol absorption and scattering phase function. The ap- 
parently low uncertainty over ocean results from the use of the wide spectral 
range of MODIS to derive the aerosol size distribution. 

A number of papers describe the validation of MODIS aerosol products, 
primarily using the operational Aerosol Robotic Network (AERONET) sunpho- 
tometers [6]. Over land, Chu et al. [3] confirm the visible x a accuracy of Ax a = 
±0.20t a ± 0.05 and a reasonable spectral correlation. Over ocean in non-dusty re- 
gimes, Remer et al. [23] confirm that MODIS retrieves x a with an accuracy Ax a = 
±0.03x a ± 0.05, and particle effective radius within an accuracy of ±0.10 fim. In an 
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oceanic dusty regime. Levy et al. [12] show similar r a errors in the mid- visible 
wavelengths, but less accuracy in determining particle radius, presumably due to 
nonsphericity effects that are not yet accounted for in the lookup tables. 

For Aqua-MODIS, the operational algorithm is the same as it is for Terra- 
MODIS. For both Terra and Aqua, improvements have been introduced for ad- 
dressing a few key issues. Over land, the reliance on dark targets has been re- 
laxed so that aerosol retrievals will be possible over a larger land surface area 
having a reflectance at 2.13 jum up to 0.4, as described earlier. This became pos- 
sible due to analysis of aircraft data on the spectral properties of land surfaces 
from desert and non-desert regions [7] and successful application to MODIS data 
[Eric Vermote, personal communication]. Over ocean, an improved cirrus detec- 
tion technique [5] has been implemented, which combines cold cloud masking, 
spatial coherence tests [17], and a new high cloud mask test based on the 1.38 (im 
band. 


VI. Total Precipitable Water Product 
MODIS has five near-infrared bands located within and around the 0.94 jum 
water vapor band for remote sensing of column water vapor amounts over clear 
land areas and over oceanic areas with sunglint. The retrieval algorithm [8], [10] 
relies on observations of water vapor absorption of near-infrared solar radiation 
reflected by the bottom surface. The algorithm uses ratios of water vapor ab- 
sorbing bands (within the 0.94 jum water vapor band) with atmospheric window 
bands at 0.86 and 1.24 fim. The ratios largely remove the effects of variation of 
surface reflectance with wavelength for most land surfaces and result in the at- 
mospheric water vapor transmittances. The column water vapor amount is de- 
rived from the transmittance using a table lookup procedure and applied to 
every 1 km pixel. The lookup tables were pre-calculated using a line-by-line at- 
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mospheric transmittance code and the HITRAN2000 spectroscopic database [24]. 
Water vapor values using the near-infrared algorithm can be determined with 
errors of about 7% based on comparisons with water vapor measurements using 
ground-based microwave radiometers. Systematic errors in the derived water 
vapor values can be further reduced if the water vapor line parameters are im- 
proved in the future. 

Fig. 7 shows a multispectral image of one Terra-MODIS granule acquired 
over the Eastern United States at 1705 UTC on June 24, 2002. Panel (a) shows a 
false color composite constructed from bands at 0.645, 0.0.858, and 0.469 gm. On 
this day, there was a strong gradient between moist air in the southeastern U.S. 
and a band of dry air extending from central Texas through eastern Nebraska 
and into southern Wisconsin. An upper-level low centered over Louisiana is as- 
sociated with the moist air throughout the southeastern United States at this 
time. The sharp boundary between this warm, moist air and the narrow streak 
of dry air coincides with a strong jet at 850 hPa that is transporting dry continen- 
tal air northward. High-resolution total precipitable water retrievals from 
MODIS can be useful in anticipating the distribution of precipitation patterns. At 
0000 UTC on June 25, 2002, seven hours after the MODIS overpass, rain was re- 
ported on either side of the dry line. Some precipitation was falling throughout 
the southeastern U.S., including thunderstorms in Illinois and Arkansas; heavy 
rain was falling in northwestern Wisconsin and Minnesota. A distinct region 
without precipitation existed along the axis of the band of low TPW observed by 
MODIS. 

The near-infrared TPW (Fig. 7b) and the thermal infrared TPW (Fig. 7c) both 
depict the strong moisture gradient very well. The 5 km IR TPW tends to smooth 
the moisture gradient and is not generating retrievals in the thin and broken 
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clouds. From the green color of the false color image, it is readily apparent that 
the surface reflectance is evident throughout the thin and broken cloud regions 
of the Gulf States, enabling selected clear sky retrievals to be obtained with this 
technique using the 1 km near-infrared bands. The near-infrared and thermal 
infrared TPW retrievals tend to corroborate each other; each provides unique 
coverage with the near-infrared 1 km product providing TPW details over land 
and sunglint regions of the ocean during the day and the thermal infrared 5 km 
product providing TPW over land and ocean during both day and night under 
clear skies. The near-infrared TPW includes a quality assurance parameter that 
indicates whether the pixel is clear or cloudy, and this cloud mask is somewhat 
different than the cloud mask used to filter the clear-sky in the thermal infrared 
algorithm. This accounts for some of the differences in coverage of the near- 
infrared and thermal infrared precipitable water parameters in Fig. 7. 

VII. Cloud Product 

The cloud product (MOD06) algorithm combines infrared and visible tech- 
niques to determine the physical, radiative, and microphysical properties of 
clouds [21]. Cloud optical thickness and effective radius are derived globally 
using 6 visible and near-infrared bands at 1 km spatial resolution. Cloud top 
properties, including cloud top temperature, cloud top pressure, and effective 
emissivity, are derived using the infrared split window and longwave CO 2 ab- 
sorption bands (both day and night) at 5 km spatial resolution. Cloud thermo- 
dynamic phase is computed at 5 km resolution using a two-band algorithm that 
includes an additional thermal band at 8.55 ;Um, and also at 1 km resolution us- 
ing a different technique based on results from the cloud mask decision tree fol- 
lowed by a bispectral threshold test, shortwave infrared tests, and finally cloud 
top temperature [21]. Finally, the cloud product contains a cirrus reflectance 
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product at a visible wavelength for use in removing cirrus scattering effects from 
the land surface, and it utilizes an additional band at 1.38 gm. Hence, the cloud 
product contains many different cloud properties derived from 14 bands in total, 
and the file size is different during the night (only cloud top properties and 
thermodynamic phase at 5 km resolution) than during the day (when additional 
1 km resolution products are included), as outlined in Table I. 

Fig. 8 shows an example of the cloud optical thickness, cloud top pressure, 
and effective radius for a daytime granule of Terra-MODIS data over the West- 
ern Pacific Ocean near the Kamchatka Peninsula on August 10, 2001 at 00:25 
UTC. The true color image in (a) shows extensive cloud cover over the Sea of 
Okhotsk, including mid-level and upper-level ice clouds, whereas the Bering Sea 
to the east of the Peninsula contains extensive marine stratocumulus clouds with 
numerous ship tracks in the southeastern portion of the image. The accompa- 
nying images show the (b) cloud optical thickness, (c) cloud top pressure, and (c) 
cloud effective radius, where we have used a different color scale for water and 
ice clouds in (b) and (d). The optically thick marine stratocumulus to the east of 
the Peninsula are identified as water clouds with optical thicknesses t c (0.65 gm) 
up to ~25. The optically thick ice clouds over the Sea of Okhotsk and around the 
southern portion of the Kamchatka Peninsula have cloud optical thicknesses ap- 
proaching 40. The ship tracks are not easily identified in the cloud optical thick- 
ness image shown here, but they result in reduced effective radii in the micro- 
physical retrievals shown in (d). For this scene, the thermal infrared and deci- 
sion tree algorithms for deriving cloud thermodynamic phase, discussed by Plat- 
nick et al [21], are in quite good agreement. Finally, the cloud top pressure for 
the cloud-filled pixels, shown in (c), clearly show that the water clouds lie pre- 
dominantly between 700 and 850 hPa, whereas the optically thick ice clouds over 
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the Sea of Okhotsk lie at altitudes above the 500 hPa level. 

VIII. Level-3 Atmosphere Products 

There are three Level-3 atmosphere products (MOD08), each covering a dif- 
ferent time interval (daily, 8-day, and monthly), that contain statistics derived 
from four Level-2 atmosphere products: aerosol (MOD04), precipitable water 
(MOD05), cloud (MOD06), and atmospheric profiles (MOD07). Statistics are 
sorted into 1° x 1° cells on an equal-angle global grid (180 x 360 cells). 

The daily product contains roughly 600 statistical datasets that are derived 
from approximately 80 scientific parameters from the four Level-2 products 
noted above. The 8-day and monthly products, identical in format, are somewhat 
larger in size than the daily product, comprising roughly 800 statistical datasets. 

For the daily product, any Level-2 granule that overlaps any part of the data 
day (0000 to 2400 UTC) is included in the computation of the statistics. There- 
fore, a particular Level-2 granule may be included in two consecutive MOD08 
daily products. The 8-day product is computed by manipulating and summa- 
rizing the daily product over eight consecutive days (half the 16-day repeat time 
of the Terra and Aqua orbits). The running 8-day intervals are reset at the be- 
ginning of each year (January 1), in concert with comparable Level-3 products 
produced by the oceans and land discipline groups of MODIS. The monthly 
product is computed by manipulating and summarizing the daily product over a 
calendar month. 

There is no separation of MODIS science parameters by orbit (ascend- 
ing/descending) node; however there is a day /night separation and/or proc- 
ess/no process decision for a number of parameters at both Level-2 and Level-3. 
All of the aerosol (MOD04) and near-infrared water vapor (MOD05) parameters, 
as well as the cirrus detection and cloud optical parameters (both a subset of 
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MOD06) are only processed for daytime (sunlit) granules. Cloud top properties 
(a subset of MOD06) are processed for both day and night; and are separated in 
Level-3 by daytime only, nighttime only, or combined day and night. Atmos- 
phere profiles (MOD07) are also processed for both day and night; but are not 
separated in the Level-3 product. 

The set of statistical summaries computed for each parameter depends on 
the parameter being considered and might include: (i) simple statistics (mean, 
minimum, maximum, standard deviation), (ii) parameters of normal and log- 
normal distributions, (iii) fraction of pixels that satisfy some condition (e.g., 
cloudy or clear), (iv) histograms of the distribution within each grid box, (v) his- 
tograms of the confidence placed in each measurement, and (vi) histograms 
and/or regressions derived from comparing one science parameter to another. 
Level-3 statistics may be computed for a subset of Level-2 pixels that satisfy some 
condition (e.g., water clouds or ice clouds). Finally, Level-3 statistics are com- 
puted by sub-sampling the Level-2 data at the resolution of the internal Level-2 
(input file) geolocation. For example, since the geolocation internal to the 
MOD06 cloud product is 5 km, every MOD08 statistic derived from MOD06 (in- 
cluding the Quality Assurance parameters) is computed from a set of Level-2 
pixels sub-sampled at that resolution (about 480 pixels instead of the 12,000 1 km 
pixels that actually fall within the 1° x 1° cell). 

In addition, histograms, joint histograms, and/or regressions derived from 
comparing one science parameter to another, are often computed for a subset of 
observations that satisfy some condition (e.g., liquid water clouds only, ice 
clouds only, etc.). The thirteen joint histograms defined in the Level-3 atmos- 
phere product (MOD08) are summarized in Table II. All joint histograms are de- 
rived from the cloud product (MOD06), ten of which are aggregated (separated) 
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by cloud phase (liquid water or ice). This group of joint histograms in computed 
from the parameters (i) cloud optical thickness, (ii) cloud effective particle radius, 
(iii) cloud top temperature, and (iv) cloud effective emissivity. In MOD08 there 
are five joint histograms produced for liquid water clouds and five for ice clouds, 
all during the daytime. The numerical boundaries used to bin the data vary by 
cloud retrieval phase, as summarized in Table II. The remaining three joint his- 
tograms are aggregated (separated) by day / night mode (daytime, nighttime, and 
combined day and night), and are computed from the parameters (v) cloud top 
pressure and (vi) cloud effective emissivity (cf. Table II). 

Although for a few parameters the set of predetermined statistics changed 
once (or perhaps twice) during the first 2 1 /2 years of MODIS data collection 
[Collection 001 (beta) and Collection 003 (provisional)], the intent was to create 
an unbroken record of fixed (unchanging) Level-3 statistical summaries, unique 
to each science parameter. For Collection 004, the first completely validated set 
of MODIS science data, the set of statistical summaries created for each science 
parameter will be fixed and unchanging. MODIS Collection 004 processing is 
slated to begin on October 29, 2002, where both a forward stream, as well as a re- 
processing stream going back to year 2000, will be produced. This will eventu- 
ally create a multiyear validated set of MODIS science data and a fixed set of 
global statistics. 

Fig. 9 shows monthly mean values of (a) cloud optical thickness, (b) aerosol 
optical thickness, and (c) precipitable water derived from the near-infrared algo- 
rithm over land and sunglint, for August 2001. The cloud optical thickness im- 
age is a combination of liquid water and ice cloud retrievals, and clearly shows 
greater optical thickness along the intertropical convergence zone as well as at 
high latitudes in the northern and southern hemispheres. Fig. 9b shows the 
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monthly mean aerosol optical thickness over both land and ocean, except for lo- 
cations where the surface is too bright to be able to retrieve the aerosol loading 
(e.g., Sahara, Saudi Arabian, and Kalahari deserts, Tibetan plateau). During 
August the biomass burning in Angola, Zambia, and the Congo are quite strik- 
ing, as are industrial emissions off China and the Indian subcontinent. In 2001 
the burning season in South America was late in starting, so August is not as dis- 
tinctive as later in September when the biomass burning season was in full 
swing. The precipitable water retrievals in Fig. 9c are derived over land and 
sunglint regions of the cloud-free ocean during the daytime using the near- 
infrared algorithm. During August there was no sunglint over the oceans for the 
Terra orbit at latitudes south of 23°S, hence no water vapor retrievals over the 
southern ocean. The large water vapor concentration in the tropics and humid 
sections of the southeastern United States and the Indian subcontinent are read- 
ily apparent. This figure compares favorably with the thermal infrared-derived 
precipitable water (not shown) that is derived over all cloud-free pixels during 
both day and night. 

These displays are the equal-area (Hammer-Aitoff) projection, which are 
available for viewing on the MODIS atmosphere web site (modis-atmos. 
gsfc.nasa.gov) along with other time periods, parameters, and statistics of inter- 
est. It should be noted that the web site contains Level-3 images on both the na- 
tive equal-angle grid as well as an equal-area grid, even though the Level-3 
product contains global statistical summaries on only the equal-angle grid. 

In addition to mean values of various parameters, the Level-3 product con- 
tains additional statistical properties including joint histograms of selected vari- 
ables, as discussed above. Fig. 10 shows a joint histogram of cloud optical thick- 
ness and effective radius for liquid water clouds over ocean for a region off the 
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Chile and Peru coast bounded by 12°-24°S and 68°-80°W for the month of August 
2001. Histograms such as this are contained within the Level-3 file for all 1° x 1° 
cells globally, and are aggregated for water clouds and ice clouds separately. 
Note the dominance of effective radii in the range of 9-14 (im, which are reason- 
able values for the boundary layer marine stratocumulus clouds prevalent in this 
region. The bin boundaries of the aggregated histogram are also shown in this 
figure. 


IX. Summary and Conclusions 

In this paper we have described each of the MODIS atmosphere data prod- 
ucts, including the frequency and size of each file, and have given examples of 
each of these data products based on observations from the Terra and, in some 
cases. Aqua spacecraft. Aqua carries a similar MODIS instrument to Terra's, but 
in a 1:30 pm ascending sun-synchronous polar orbit, rather than the 10:30 am de- 
scending orbit of Terra. The same data products are produced from Aqua as 
from Terra, thus enabling studies of diurnal variability of cloud, aerosol, and wa- 
ter vapor properties. Due to differences in the performance of individual detec- 
tors, cross talk, and gain settings, there are some subtle and not so subtle differ- 
ences in several of the data products. Of particular significance is the known 
failure of many of the individual detectors in the 1.64 gm spectral band that 
makes it impossible to retrieve cloud effective radius using this band. Otherwise, 
preliminary examination of the Aqua products shows a similar accuracy to the 
Terra data products. Therefore, the experience in the precise characterization of 
the MODIS spectral bands on Terra that was achieved after the first year of op- 
eration yielded similar precision on the Aqua spacecraft after only 3 months in 
orbit. The location of sunglint is considerably different in the afternoon orbit, so 
the areas of the ocean that are impossible to retrieve water vapor concentration 
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from using the near-infrared bands (like the southern oceans during August 
2001, as shown in Fig. 9c) will now move to the northern hemisphere. Hence, be- 
tween Terra and Aqua it should be possible to obtain high spatial resolution pre- 
cipitable water vapor concentrations globally. 

All data produced by the algorithms described in this paper are available 
from the Goddard Distributed Active Archive Center, and may be ordered via 
the EOSDIS Data Gateway (redhook.gsfc.nasa.gov/~imswww/pub/imswel- 
come) or Terra WHOM (daac.gsfc.nasa.gov/data). Quicklook images of all 
Level-IB true color and Level-3 data products are available from the MODIS at- 
mosphere web site at modis-atmos.gsfc.nasa.gov. In addition, this site contains 
(i) software tools for locating granules, doing spatial subsetting, and visualizing 
various scientific data sets, (ii) extensive references and descriptions of the algo- 
rithms used to process the data, and (iii) calendars of data availability. 
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Table I 

MODIS Atmosphere Data Products, Responsible Investigator, Processing 
Level, Product ID, Spatial Resolution, File Frequency, and File Size. 


Data Set 

Responsible 

Investigators 

Process- 
ing Level 

Product 

ID a 

Spatial 

Resolution 

(km) 

File 

Frequency 

File Size 

Pixel-level products 
Qoud mask 

Ackerman, 

Menzel 

2 

MOD35 

0.25, 1 

288 /day 

47.4 MB 

Aerosol product 

Kaufman, 
Tanrd, Remer 

2 

MOD04 

10 

139/ day 

12.0 MB 

Total precipitable water 

Gao, Kaufman, 
Tanr£, Menzel 

2 

MOD05 

1 (NIR), 
5 (TIR) 

288 /day 

20.2 MB (day) 
3.6 MB (night) 

Qoud product 

King, Platnick, 
Menzel, Gao 

2 

MOD06 

1/ 5 b 

288 /day 

69.6 MB (day) 
14.1 MB (night) 

Atmospheric profiles 

Menzel 

2 

MOD07 

5 

288 /day 

32.2 MB 

Gridded time-averaged products 
Level-3 atmosphere prod- 
ucts 

King, Hubanks, 
Pincus 

3 

MOD08 

l°xl° 

1/day 
1/ 8-day 
1/ month 

440.9 MB 
810.6 MB 
810.6 MB 


a MODxx products apply to Terra and MYDxx products apply to Aqua. 
b Daytime retrievals of cloud optical and microphysical properties are performed at 1 km 
resolution, while daytime and nighttime cloud top altitude and cloud top pressure are 
processed at 5 km [21]. 
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Table II 

MODIS Cloud Parameters that Appear in Joint Histograms. 


Parameter 1 

Day/ 

Night 

Bin Boundaries 

Parameter 2 

Bin Boundaries 

Liquid water clouds 





Optical thickness 

Day 

0, 2, 4, 6, 8, 10, 15, 
20, 30, 40, 50, 100 

Effective radius (jum) 

2, 4, 6, 8, 10, 12.5, 15, 17.5, 
20, 25, 30 

Optical thickness 

Day 

0, 2, 4, 6, 8, 10, 15, 
20, 30, 40, 50, 100 

Qoud top tempera- 
ture (K) 

230, 235, 240, 245, 250, 255, 
260, 265, 270, 275, 280, 285, 
290 

Optical thickness 

Day 

0, 2, 4, 6, 8, 10, 15, 
20, 30, 40, 50, 100 
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FIGURE LEGENDS 

Fig. 1. MODIS atmosphere data processing architecture and products (adapted 
from [20]). 

Fig. 2. Cloud mask over Western Australia on June 5, 2001. Panel (a) is a true 
color composite of one MODIS granule using bands at 0.645, 0.555, and 
0.469 jum. Panel (b) is the cloud mask derived for this scene, showing 
the various confidence levels of cloud identified in this scene. 

Fig. 3. Atmospheric temperatures retrieved from Terra-MODIS radiances on 
May 22, 2001 at (a) 850 hPa, (b) 700 hPa, and (c) 500 hPa. 

Fig. 4. Cross-section of (a) temperature and (b) mixing ratio along a cross sec- 
tion in the southeastern United States (white line in panel (c)), which 
shows the thermal infrared-derived total precipitable water vapor for an 
Aqua-MODIS granule on September 4. 2002. (d) The associated GOES-8 
Sounder depiction of the total precipitable water vapor. The white bars 
in (a) and (b) correspond to regions obscured by cloud. 

Fig. 5. Comparison of total precipitable water from MODIS thermal-infrared 
regression (solid circles), GOES-8 (open circles), and radiosonde (open 
diamonds) with the SGP ARM-CART microwave radiometer (MWR) in 
millimeters. 80 MODIS cases from April 1 to September 1, 2002 are 
shown in the comparison. 

Fig. 6. Aerosol optical thickness over the Eastern United States on May 4, 2001. 
Panel (a) is a true color composite of one MODIS granule, showing cloud 
cover over northern New York, central Florida, and the nearby Atlantic 
Ocean. Panels (b) and (d) show the aerosol optical thickness at 0.56 jum 
and the effective radius derived for this scene, showing pollution along 
the Ohio Valley and transport to the north Atlantic. Panel (c) shows the 
ratio of aerosol optical thickness of the fine to the total aerosol optical 
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thickness. Gray areas denote areas where no retrieval is possible due to 
sunglint, bright land, and clouds. 

Fig. 7. Water vapor over the eastern United States on June 24, 2002. Panel (a) is 
a false color composite of one Terra-MODIS granule using bands at 
0.645, 0.858, and 0.469 jum, showing cloud cover over the upper Midwest 
and Mid-Atlantic regions of the United States. Panel (b) is total precipi- 
table water derived from the cloud-free pixels over land using the near- 
infrared algorithm. Panel (c) is the total precipitable water derived from 
the thermal infrared algorithm. 

Fig. 8. Cloud properties over the western Pacific Ocean off the Kamchatka 
Peninsula on August 10, 2001. Panel (a) is a true color composite of one 
MODIS granule, showing marine stratocumulus clouds with ship tracks 
as well as upper level ice clouds. Panels (b) and (d) show the cloud opti- 
cal thickness and effective radius derived from all cloudy pixels, where 
we have used a separate color bar to denote clouds processed as ice and 
water clouds. Panel (c) shows the cloud top pressure for all clouds in 
this scene. 

Fig. 9. Level-3 monthly averages of (a) cloud optical thickness, (b) aerosol opti- 
cal thickness, and (c) near-infrared precipitable water for August 2001. 
These monthly products are produced at a 1° x 1° latitude-longitude grid 
worldwide. 

Fig. 10. Histogram of the joint probability of cloud optical thickness and effec- 
tive radius derived from the Level-3 joint atmosphere product for 
August 2001 for a region off the coast of Peru and Chile bounded by 12°- 
24°S and 68°-80°W. All liquid water cloud retrievals over ocean in this 
latitude-longitude box have been aggregated. The grid pattern reflects 
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the bin boundaries in this Level-3 data product. 
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Fig. 1. MODIS atmosphere data processing architecture and products (adapted from [20]). 
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Fig. 2. Cloud mask over Western Australia on June 5, 2001. Panel (a) is a true color composite of one 
MODIS granule using bands at 0.645, 0.555, and 0.469 #«n. Panel (b) is the cloud mask 
derived for this scene, showing the various confidence levels of cloud identified in this scene. 
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Fig. 3. Atmospheric temperatures retrieved from Terra-MODIS radiances on May 22, 2001 at (a) 
850 hPa, (b) 700 hPa, and (c) 500 hPa. 
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Fig. 4. Cross-section of (a) temperature and (b) mixing ratio along a cross section in the southeastern United States 
(white line in panel (c)), which shows the thermal infrared-derived total precipitable water vapor for an 
Aqua-MODIS granule on September 4. 2002. (d) The associated GOES-8 Sounder depiction of the total 
precipitable water vapor. The white bars in (a) and (b) correspond to regions obscured by cloud. 
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Fig. 5. Comparison of total precipitable water from MODIS regression (solid circles), GOES-8 
(open circles), and radiosonde (open diamonds) with the SGP ARM-CART microwave 
radiometer (MWR) in millimeters. 80 MODIS cases from April 1 to September 1, 2002 are 
shown in the comparison. 
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Fig. 6. Aerosol optical thickness over the Eastern United States on May 4, 2001. Panel (a) is a true 
color composite of one MODIS granule, showing cloud cover over northern New York, 
central Florida, and the nearby Atlantic Ocean. Panels (b) and (d) show the aerosol optical 
thickness at 0.56 fun and the effective radius derived for this scene, showing pollution 
along the Ohio Valley and transport to the north Atlantic. Panel (c) shows the ratio of 
aerosol optical thickness between the fine and coarse particle mode. Gray areas denote 
areas where no retrieval is possible due to sunglint, bright land, and clouds. 




KING et al.: ATMOSPHERIC PROPERTIES FROM MODIS 


41 



a) R(0.645, 0.858, 0.469) 


b) Near-Infrared Precipitable Water 



0 12 3 4 5 6 

Water Vapor (cm) 


c) Thermal Infrared Precipitable Water 



Fig. 7. Water vapor over the eastern United States on June 24, 2002. Panel (a) is a false color 
composite of one Terra-MODIS granule using bands at 0.645, 0.858, and 0.469 ftm, showing 
cloud cover over the upper Midwest and Mid-Atlantic regions of the United States. Panel (d) 
is total precipitable water derived from the cloud-free pixels over land using the near- 
infrared algorithm. Panel (c) is the total precipitable water derived from the thermal infrared 
algorithm. 
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Fig. 8. Cloud properties over the western Pacific Ocean off the Kamchatka Peninsula on August 10, 
2001. Panel (a) is a true color composite of one MODIS granule, showing marine 
stratocumulus clouds with ship tracks as well as upper level ice clouds. Panels (b) and (d) 
show the cloud optical thickness and effective radius derived from all cloudy pixels, where 
we have used a separate color bar to denote clouds processed as ice and water clouds. Panel 
(c) shows the cloud top pressure for all clouds in this scene. 
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Fig. 9. Level-3 monthly averages of (a) cloud optical thickness, (b) aerosol optical thickness, and (c) 
near-infrared precipitable water for August 2001. These monthly products are produced at a 
1° x 1° latitude-longitude grid worldwide. 
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Fig. 10. Histogram of the joint probability of cloud optical thickness and effective radius derived from the Level-3 joint 
atmosphere product for August 2001 for a region off the coast of Peru and Chile bounded by 12°-24°S and 68°-80°W. 
All liquid water cloud retrievals over ocean in this latitude-longitude box have been aggregated. The grid pattern 
reflects the bin boundaries in this Level-3 data product. 


